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Abstract 

The characteristic functions of multivariate Feller processes with 
generator of affine type, and with smooth symbol functions have an 
explicit representation in terms of powers series with rational number 
coefficients and with monoms consisting of the symbol functions and 
formal derivatives of the symbol functions. The power series repre- 
sentations are convergent globally in time and on bounded domains of 
arbitratry size. Generalized symbol functions can be derived leading to 
power series expansions which are convergent on unbounded domains. 
The rational number coefficients can be efficiently computed by an in- 
teger recursion. As a numerical consequence characteristic functions of 
multivariate affine processes may be computed directly from the sym- 
bol function avoiding the generalized Riccati equation (an observation 
first made in [3] in a more general context). 

1 Introduction 

Bochner told us to investigate Markov processes by the symbols of their 
associated infinitesimal generators. Recall the simple connection. From the 
probabilistic point of view we may start with a regular Feller process 

with values in M. d . Then the function 

E x ( e i ( x ?- x )<\ - l 

is called the symbol of the process. The function 

ct(x,0=E x (e i{X ?- xH ) (2) 
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entering the definition in ([I]) is called the characteristic function. Now, if 
(Tt)t>o is called the semigroup associated with ((Xf) t >o, P x ) xGR „, then we 
have 

T t f(x) = E x (/(Xf)) = (27T)-"/ 2 / e^c t (x,Of(m, (3) 

JR n 

where ? denotes the Fourier transform. It follows that the generator takes 
the form 

Au(x) = lim TMX) ~ U{X) = -(2vr)^/ 2 / e**a(z, (4) 

Hence the symbol completely characterizes the process. However, recovering 
the characteristic function of a process from its symbol is far from trivial. 
Especially it is far from trivial to say under which condition the symbol 
of a given partial integro-differential operator corresponds to a Markov pro- 
cess. Hormanders standard class of symbols of pseudo differential operators, 
namely functions (x, £) — > a(x, £) which satisfy 

ia^fa(x,e)i< Ca/3 yrw m " p|Q|+5|/31 , ( 5) 

where 0<5<p<ldo not fit very well in order to analyze this ques- 
tion. Symbolic calculus was generalized in pQ, and further generalized by 
[6], where the function yl + |£| 2 is replaced by more generalized weight 
functions. In this spirit weight functions were constructed such that the 
corresponding symbol class corresponds to a class of Feller processes (for 
an overview of more recent results cf. [7J and references therein). Well it is 
clear, that symbolic calculus cannot fully characterize the class of Feller pro- 
cesses at least in its usual form because even in the class of Levy processes 
you find symbols which are not differentiable but correspond to fractional 
Laplacians for example. (Well, some symbolic calculus for non-differentiable 
symbol functions has been developed, but seems not suitable in the present 
case). However, the class of Feller processes with differentiable symbols is 
quite considerable and for some fields of applications, especially finance, it 
seems to be satisfying. Well within this class the next question to be con- 
sidered is how we can construct a certain class of Feller processes from their 
symbols, and the relation outlined above leads us to aim at the characteristic 
function of the process. For processes X where for each t > Xt has an 
infinitely divisible distribution (such processes correspond Levy processes in 
law), an explicit formula for the characteristic function is known in form of 
a Levy-Kintchine formula, i.e. 

E \ e i{u,Xf-x)\ ._ e -\(Au,u)t+i(r,u)t+t ^ d {e^)-l-i(u,x)l D {x)u{dx)) ^ ^ 

where A is a symmetric nonnegative dxd matrix, r € and v is a measure 
on W 1 satisfying 

^({0}) = 0, and / min { \y\ 2 , l} v(dy) < 00. (7) 
JR d \{0} 
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For affine processes explicit formulas for the characteristic function are 
known only in special cases (cf. the example in [3] for the case of affine 
processes without jumps). In the general case the computation of the (con- 
ditional) Fourier transform of an affine process is shown to be reducible to 
solving systems of general Riccati differential equations (cf. [5]). Analyti- 
cal solutions of the generalized Riccati differential equation systems are not 
known, and even numerically they are often difficult to solve, especially due 
to quadratic terms involved. In [3] it is shown that for a class of affine equa- 
tions a basis of analytic vectors can be found. In particular, this leads to 
recursive representations of characteristic functions in terms of constituents 
of the generator avoiding the generalized Riccati equations. In a special 
univariate case an explicit formula is derived. 

In this paper we concentrate on the characteristic function and show 
that an explicit series in terms of the symbol and formal derivatives of the 
symbol with rational number coefficients can be derived. Convergence of 
the representation is shown. 



2 Explicit formulas for characteristic functions of 
affine processes 

For affine process X°' x with X°> x = x we search for an explicit representation 

of p(t,x,u) = E[e mXt ' ). Here we denote e tux := e^i=^ u i x i for simplify of 
notation. We assume that the function p is a global solution of the Cauchy 
problem 

%{t,x,u) = Ap(t,x,u), 



(8) 



p(0,x,u) = exp(mx), t > 0, x € fl C 
where is some domain and 

d 2 

"x) 



mix) = ^£m*)^/(*) + £m*)^;/(3 

ij 1 ] 



+ 

/R d \{0} 



df 

f(x + z) - f(x) - -q^{x) ■ zl D (z) 



v(x, dz). 



Here we assume that the equation is semielliptic, i.e. 

(aij) > 0, 

and we assume that the operator is closed on some appropriate function 
space of its domain- in the case of a bounded Q this is not a problem, but 
in general additional conditions on the measure v may be needed (cf. [3] 
for a discussion of that point). Furthermore 1^ is a indicator function of 
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a bounded set D, aij(x), bi{x) are affine functions, and v(x,dz) is a Borel 
measure which is affine in x, i.e. we have the representations 

aij(x) = dijo + Ya=i a iji x i 

bi(x) = b i0 + Ya=i baxi (9) 

u(x,dz) = ln o u (dz) + Yli=i x i 1 o, l ^i{dz), 
where fioi^! C I'', 1 < I < d. We define the symbol function 



a(x,iu) := exp (— iux)A[expiux] = — -^2aj k (x)ujU k + bj(x)iuj+ 

jk j 

exp (— iux) I [exp(zii(x + z)) — exp{iux) — iuex.p(iux) • zljj(z)] v(x, dz) 



ajk(x)ujUk + bj(x)iuj + / [exp(inz) — 1 — iu ■ zlr>(z)] u(x, dz) 

2 # T jRd \w 

Since a is affine with respect to x it is useful to define 

aiiiu) := d Xl a(x, iu) = -\ J2 jk a jk iuju k + J2j b ji iu j 

(10) 

+ Ef=i Iri\{ 0} [exp(raz) - 1 - iu ■ zl D (z)} vi(dz) 
In the case d = 1 we define 

a\{iu) := §^<j{iu). (11) 

Remark 2.1. In general we denote multivariate derivatives of order a = 
(ax, ■ ■ ■ , ad) of the symbol function with respect to the symbol variable £ by 
d^a. Note that for \a\ > 3 multivariate derivatives of the symbol function 
are bounded, i.e. \d^a\ < oo if and only if | J R d\/ ]. z a u(x, dz)\ < oo. 

As a main result of the present paper we derive power series formu- 
las for the characteristic functions of affine processes with rational number 
coefficient and (multivariate) monomials of the form 

y k / \ H k 

- (iu) 



U\ ejl < k Uf =1 [d^a(x,iu)) 

where the index j denotes an enumeration of all multiindices e, of dimension 
d and of order k where k varies over all nonnegative integer, i.e. e,- = 
(ejx, ■ ■ ■ ,ejd) and X^f=i e ji = ^, and a| and /3j denote nonnegative integers 
enumerated in accordance with the multiindices ej. Obviously, we have 
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Proposition 2.2. The symbol function 

a : n x iR d -> C, SlCl d open, 

a(x,iu) = -l^2 jk a jk (x)uju k + £\ bj(x)iuj (12) 

+ / R(J [exp(mz) — 1 — iu ■ z1d(z)] v(x,dz) 

is globally bounded in the spatial variables x if and only if either 0, is 
bounded or bji,aj k i = for 1 < I < d and the measure v = ln vo(dz) + 
Ylt=i x i^Qi u i(dz) has bounded Qi, 1 < I < d. 

First we shall derive formulas which are valid locally. Then in a second 
step we globalize the result to bounded spatial domain and global time (in 
the spirit of [3]). More precisely, in order to obtain a representation which 
is globally in time we make use the following transformation 

t : [0,1) -> [0,oo), 

(13) 

t(r) :=/51n (tan(f + f r) . 

Note that in the transformed time coordinates we get an equivalent Cauchy 
problem, namely 

%^(t, x, u) = f s . // +2 r r \ (t, x, u), 

2 (14) 
p*(0, x, u) = exp(iux), t € [0, 1) ifSlC M d . 

Remark 2.3. In order to keep notation simple we shall abbreviate 

4 szn (| + |rj 

Remark 2.4. In [3] The transformation t(r) = — /31n(l — r) is used for similar 
purposes. Our choice is guided by computational/numerical considerations 
which will be pointed out in a subsequent paper. 

It makes sense to consider the univariate case first, because even that 
case is involved. 

2.1 The univariate case 

We introduce the following multiindex notation. For each positive integer k 
we denote 

a k = (a k ,--- ,ati), P k = ■ ■ ■ , tt-i)- (16) 
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where the entries ct k ,f3 k are nonnegative integers. For any positive integer 
k the projection of a /c-tuples a k (resp. f3 k on a Z-tuple for an nonnegative 
integer I < k will be denoted by IT^ with 

U k a k = (<& • • • ,«ti), (resp. Hf/3 fc = (/3 fc , • • • (17) 
Adding or subtracting integers to such multiindices is denoted as follows: 



a k + L 



(a 



0) 



i Q! j -\- ij, • 



> a fc-l 



)■ 



(18) 



Similar for f3 k and for subtraction, i.e. we have 



a k - 



"j - K,-- - ,ctj-ij,--- (19) 

In our expansion only tuples with nonnegative entries will appear. However, 
we do not exclude negative integer entries from the beginning in order to 
keep the rules simple. Certain constants depending on the fe-tuples a k ,fi k 
will just be defined zero if negative entries occur, and then they do not 
contribute to our expansion. We have the following result 

Theorem 2.5. Locally, the following representation holds: 

p{t,x,u) 



exp{mx) x 



k-l 



1 + Y^k>l,(a k ,/3 k )€M k C (a fe ,/3 fe ) n j=0 



dia(x,iu)) 3 (diai(iu) 



3 +k 



(20) 



where 



M, 



(a k ,(3 k ) 



k-l 

£ 

3=0 



k-l 



a k + p k )=kkJ2^> 



3=0 



k-l 



(21) 



Here, fork>2 c (afc)/3fc ) = if {a k ,(3 k ) £ M k , and 
c (i,o) = 1) c (o,i) = 
and for k > 2 

1 



C( a fc+i )/3 fc+i) (^pjyy 
C (a k + 1 ,/3 k + 1 ) = (jfe+l)! 

Ek—1 



if a k+1 



(22) 



/n fc o fe+1 -ij+A fc 



) c (n k a 



k + 1 -i j +x k ,n k /3 k + 1 -x k ) 



Here, the natural restriction holds that X k = (Aq, • • • , A^_ 1 ) are k-tuples such 
that 

Il k a k+1 - lj + A fc € M k . (23) 

(Alternatively we could define cr a h^k-\ with pairs of k-tuples (a k ,f3 k ) not in 
M k to be zero, of course). 
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Remark 2.6. The recursion is very easy to compute. Let us compute the 
approximation for the first few terms in full expansion (up to order 3 of (|22p). 
Let dk denote the sum of terms of order k in time. We have 

d Q = 1 

d\ = cr(x, iu) 

d2 = \a{x,iu) 2 + \{d^a{x,iu))ai{iu) (24) 

d>3 = ^a(x,iu) 3 + ^a(x, iu)(d^a(x, iu))ai(iu) 

^(d^a(x,iu))(d^(ai(iu))ai(iu) + ^ u?|cr(x, iu)J o\{iu) 2 . 

Well, higher order terms are easy to compute as well. The formular is coded 
by an affine triangle of rational numbers with fcth row all numbers cr a k t pk\ 
where a k , f3 k are the fc-tuples. We call that row the univariate affine triangle 
of the caracteristic function. The first row (terms of order 1) contains just 
one number 

C(i,o) = 1 (25) 
(note that the tuple (0, 1) £ Mi). The second row contains two numbers 

1 1 

C(( 2) o),(o,o)) = 2> c ((o,i),(i,o)) = 2 ( 26 ) 

The third row contains four numbers 

c ((3,0,0), (0,0,0)) = 6 c ((l, 1,0), (1,0,0)) = 2 

(27) 

c ((0,l,0), (1,1,0)) = 6 c ((0,0,2), (2,0,0)) = 6' 

and so on. The growth of the number of numbers in the rows of the affine 
triangle of the characteristic function is exponentially. However, with a 
simple computer program computation of order 20 (about a million numbers) 
is no problem leading to formulas of unprecedented precision. 

Remark 2.7. Note that in the second sum of[22]the cardinality of the second 
sum is 

CI!; 1 ) 

a information which is useful for consideration of convergence. 

Remark 2.8. The representation is different from the special case considered 
in [3] . In that case there is only one affine function in front of the integral 
operator, and we used this special structure to shift the affine functions 
within the symbol expansion. This is not possible in the general univariate 
affine case considered here. 
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Next we 'globalize' the result (quite in the spirit of [3] by considering 
transformed variable). Note that in the case of transformed variables t{r) = 
— /31n(tan(| + jt)) we have symbol functions which depend on time, i.e. 

Op : [0, T T )x!lx iR d -> C, 



a p (T,x,iu) := p(r)a(x,iu) 

However, we let the symbol variable £ vary on the whole domain iM. d 
for the class of bounded symbol functions with bounded derivatives. As we 
observed, this is the case for processes with bounded moments and constant 
diffusion and drift coefficients. 

Theorem 2.9. Let f2 C K. be some bounded domain of arbitraty size and 
let T G (0, oo) some horizon. First assume that the symbol function a is 
bounded on fi x iM. d . Then there is a a parameter value /3 > such that on 
[0, T T ] x x M. d the following representation holds: 

P P (r, x, u) = exp(iux) (l + Ylk>i d U T > x > ^)r fc ) , (29) 

where 

d* k {r,x,iu) : = 

(30) 

EjLl T,(m,l)eM- kj n r=l (P mr ) lr d EU mr ^ ^ Ei=S.+V)! ' 

M T m := \(m,l) := (mi,h,--- , mj, lj) G J](mr + U-) = k 1 , (31) 



r=l 



and for I = X^r=i m r 



di(x,iu) = ^2 c (a fc ,/3 fe ) n i=o [&a{x,iu)\ 3 [d^ai(iu)J (32) 

(o fe ,/3 fc )6M fe 

with the same recursion for Cf a k^k\ as before. Furthermore (p m )^ denotes 
the Ith derivative of the mth power of the function p. Explicit formulas for 
that are easily available by elementary calculation. 

Note that this implies a representation of the characteristic function 
which is global in time by replacing r by r{t) in the representation above. 

In the general case we have 

Corollary 2.10. Let fi, 5 C R be bounded domains of arbitraty size and 
let T G (0, oo) some horizon. Then there is a a parameter value (3 > such 
that on [0, T T ] x £1 x S the representation described in in 129\) . L30\) . 13 
MB holds. 
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The latter corollary shows that univariate affine processes with smooth 
symbols can be approximated on finite domains with accuracy of any prac- 
tical interest by the power series described in theorem 2.9.. However, in 
Section 4 we shall define generalized symbol functions and describe general- 
izations of theorem 2.9 where we describe power series in terms of generlized 
symbol functions and their dervatives which converge on the whole domain 
(i.e. Q C 1™ may be unbounded or equal to M. d and S = i]R d in corollary 
2.10). 

2.2 The multivariate case 

The difference between the univariate and the multivariate case is not that 
big as one may think. Well, we shall minimize this difference using a certain 
trick (which is also useful when programming the formulas of this article). 
For each positive integer k we enumerate the d-dimensional multiindices in 
a list 

ej = (e,i, • • • , e jd )\ ^ e jr = k k j € {0, • • • N k - 1} I (33) 
r=l J 

where N k is a natural number which euals the cardinality of the (i-tuples of 
order k. In addition we need some 'fc-tuples of ci-tuples' For each positive 
integer k we denote 

= (^d0> " ' j \ffc-i)) Pd = (fidoi ' ' ' ^d(k-i))- 
where the entries are (i-tuples, i.e. 

^dj = (^dj'li " > ^djd^j ) Pdj = (Pdjli ' ' " ) PdjdJ > (35) 

and with nonnegative numbers ^dji^dji ^ — ^ — d- Adding or subtracting 
integers to such multiindices is defined as before. Finally we have to adjust 
the tuples and the projection operators for tuples a bit. Instead of a k in the 
univariate case we now have tuple 

a k = (a , ■■■ , ajVfc-i)- 

For all k £ N we define the projection operator IljVj. on the N^+i tuples of 
nonnegative integers which is determined by the relation 

IIjV fe (Q!0,-- - = (OO,--- , a N k )- 

Furthermore for each k G N we define 
We have the following result 
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Theorem 2.11. Locally, the following representation holds: 
p(t,x,u) = exp(rax)(l + T,k>l,{a N k ,/3%)eM$ V^)* 



(36) 

n^Dt! {d^a(x,i U )Y j (fl£' 0| («,))"* , t*) > 

where 

N k -1 d Ar fc -l d fc-i 

M*=<j(cA/3 d fc ) E £ («" + 40 = fc & EE^>E a 

j=0 J=l j=0 1=1 j=l 

Here, fork>2 c {ctkA) = if (a k ,Pk) £ M$, and 
cnjo,— ,o)) = l( a ^ other tuples of order 1 eguaZ zero 
and /or k > 2 

1 • _£ ^*fc+l ; i 

C (a JV fe+i,^+ 1 ) = (k+W %i "0 = fc + 1 

c (a"*+i,/3* +1 ) = (fc+T)T X 

C(n^ a ^+i-i i+ A^ ) n^*+ 1 -A|)- 

Proof. Polynomial growth of the derivatives witb respect to the symbol vari- 
able means that the symbol function is analytic (this can be checked easily 
by estmating the remainder term of the multivariate Taylor formula) . Hence 
we have 

a{x,iu) = -Ef J =i%'fc^fc + Ej=iM^ + £|<5|>i ! ^M 5 

(39) 

=:£| 5 |>l^(-) 5 , 
where 5 denotes d-tuples with nonnegative intgegers, and for \5\ > 2 

m 5 (x) : = / z S u(x,dz). (40) 
Jm d \{0} 

Plugging in the ansatz p(t,x,u) = exp(iu T x) (l + dk(u, x)t k ) leads to the 
recursive formula 

4+1 = Fpr exp(-inx) J2s>o -ft® 6 (exp(mx)d fc ) (41) 



(37) 



(38) 
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Using multivariate Leibniz derivative rule we get 

4+1 = kh; exp(-iux) *£ 5 > Q ™ J2e<8 © (^) <5 " e exp(iux)d e (d k ) 

(42) 

= Y^\e\<k 9 l a ( X ' iu )?. d x( d k) 

□ 

Remark 2.12. The recursion is very easy to compute. Let us compute the 
approximation for the first few terms in full expansion (up to order 3 of (|22p ). 
Let dk denote the sum of terms of order k in time. We have 

d = 1 

d\ = a(x,iu) 

d.3 = j$cr(x, iuf + ^^2 l=1 a{x,iu)(d^a(x,iu))ai(iu) 

+gi£*,j=i (^eU^OMu)) ^W 2 - 

and so on. 

Again this theorem may be globalized leading to 

Theorem 2.13. Let Q C M. d be some bounded domain of arbitraty size and 
let T £ (0, oo) some horizon. First assume that the symbol function a is 
bounded on fi x iM. d . Then there is a a parameter value (3 > such that on 
[0, T T ] x Q, x M. n the following representation holds: 



P P (r, x, u) = exp(iux) ( 1 + Y,k>i d t( T ' x > iu ) rk ) > ( 44 ) 

where 

d* k (T,x,iu) :- 

E.Li 52(m,l)eMz UJ n r=i {p mr ) {lr) d Tj (x, iu) ^F r Tl m " 



(45) 



.3=1 ^(m,l)eM£ k . -r=l VP ; E 3 r= i"V l ' ; Ei=i("V+'r)!' 

A^Sfci : = IkO := (mi, Zi,-- - ,m j ,l j )en 2 j \J2(^r + lr) = k\, (46) 



r=l 
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and for s = J2i=i m r 

d s (x,iu) = Y,(a N s,l3 s d )£M° C (a N s,l3°) X 

with the same recursion for cr a k pk\ as before. Furthermore (p m )^' denotes 
the Ith derivative of the mth power of the function p. Explicit formulas for 
that are easily available by elementary calculation. 

Note that this implies a representation of the characteristic function 
which is global in time by replacing r by r(t) in the representation above. 

In the general case we have 

Corollary 2.14. Let Q, S C R be bounded domains of arbitraty size and 
let T € (0, oo ) some horizon. Then there is a a parameter value (3 > such 
that on [0, T T ] x S7 x S the representation described in in HPj) . [J5\ ), (4°V > 
g3> holds. 

The latter corollary shows that univariate affine processes with smooth 
symbols can be approximated on finite domains with accuracy of any prac- 
tical interest by the power series described in theorem 2.9.. However, in 
Section 4 we shall define generalized symbol functions and describe general- 
izations of theorem 2.9 where we describe power series in terms of generlized 
symbol functions and their dervatives which converge on the whole domain 
(i.e. $7 C R n may be unbounded or equal to R d and 5 = iR d in corollary 
2.10). 

3 Convergence by counting terms 

The structure of the explicit formulas makes it possible to prove convergence 
by counting the number of terms in the explicit expansion. We first have a 
look at the univariate case and then will generalize to the multivariate case. 

Let with 1 < k < n, k,n G N denote the number of terms of order 
n in time, i.e. coefficients of t n in the expansion (|20p . and of order k in the 
spatial variables, i.e. the monoms of the power series (|20p ) of order k, where 
k = ki + k2 + k%, and where k\ is the exponent of drift coefficient functions or 
derivatives of drift coefficients involved, and ki is the exponent of diffusion 
coefficient functions or derivatives of diffusion coefficient functions involved, 
and A?3 is the exponent of integral terms with an affine jump measure or 
derivatives of integral terms with an affine jump measure involved. Note 
that in difference to the expansion (|2U|) we collect all terms of order k in 
spatial variables for some time order n, i.e. we make no difference between 
monoms where mutually different derivatives of the symbol function occurr. 

12 



(47) 



Then we may determine that numbers in an integer recursion which we 
call the 'Affine Counting triangle'. In the following triangle the number of 
the row denotes the n, and the number of the column denotes the k of U r £. 
The first rows of the Affine Counting triangle are 

1 

1 1 

1 6 10 3 (48) 

1 10 34 45 4 



Elementary consideration lead us to 

Lemma 3.1. The number of terms of order n in time and of order n — k 
w.r.t. the spatial variables in the explicit formula is given by the recursion 

n™ = 1; n € N 

(49) 

K-k = Eto r k -7 l )K-U *, n € N, 1 < k < n 
Lemma 3.2. An elementary induction argument leads us to 

n-l 

Y,K-k<n\ (50) 
fc=o 

Remark 3.3. Indeed, the latter result is a rough estimate. Defining 

n-l 

R k = ^2K_ k <n\, (51) 

from the first five rows of the Affine Counting triangle we see that 
Ex - i Ez - 21 _ i R 3 _ 6 _ -, 

1! ' 2! 2! ' 3! 3! ' 

(52) 

Ri — 20 _ 5 ife _ 94 — 47 
4! — 4! — 6' 5! — 5! — 60' 

It can be shown that -§f 4 as k f oo, but I have not determined the exact 
ratio of convergence yet. Such an analysis may be important for numerical 
and computational purposes. However the rough estimate by the faculty is 
enough for our purposes. 

We have 

Lemma 3.4. The power series S20\) converges locally to the characteristic 
function of the process . 



13 



Proof. First we consider the case where the process has a bounded sym- 
bol function. Then we have less then k\ monoms of order k each one has 
the factor rr and consists of a finite product of the symbol function and 
drivatives of the symbol function. Moreover we have seen that if the multi- 



plicand yd^o~(x, iu)j 3 with < j < k — 1 in such a symbol is accompanied 

by a multipicand which is a powere of o~i{iu) with same exponent. Since 
&j;0~(x,iu) < c^a(x,iu) each momom is bounded by c k for some < c < 1. 
This means \dk\ < c k and for r € [0, 1) the whole series is absolutely bounded 



4 Generalized symbol function 

In general it is known that the characteristic function of regular affine pro- 
cesses (X, (P x )xett) with associated semigroup (T t )t£R + have a representa- 
tion of the form 



T t exp(iux) = exp (<f>(t, u) + ip(t, u)x)) (53) 

where the functions 4> an d ip are determined by generalized Riccati equations 
(cf. [5] for details). Well these are difficult to solve (even numerically), but 
in some special cases it can be done and was done. In such cases we may use 
the explicit formula and define generalized symbol functions from it. The 
the procedure is as follows. Assume that for some Cauchy problem with 
reduced symbol function you can solve the Cauchy problem For example 
you may think of some equations of type 

' f (t, x, u) = (i Zij "ui-rl-j^r + Ei Hx)ir) P (*> *> u ) = °» 

(54) 

p(0,x,u) = exp(iux), t > 0, ie!lc M d , 

where is some domain and a™ and 6, are affine function as in . Examples 
for this case are given below. However there are also examples of explicit 
solutions in form of (|53p if there are jump terms with Levy measure (cf. [?]). 
Anyway let us assume that we have an explicit solution of form (again we 
write X^f=i ^o(i) u)ixi for simplicity of notation) 

exp (4> (t, u) + ip (t, u)x) (55) 

of the Cauchy problem 

^(t,x,u) = A p(t,x,u), 

(56) 

p(0, x, u) = exp(mx), t > 0, x G Q C M. d , 
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where SI is some domain and Aq is some operator 



1 3 7. ^ 



I.I 



+ 



df, 



f(x + z) - f(x) - q^(x) ■ zl D (z) 



v°(x, dz). 



Keep in mind that most examples will be of the form (|54p or maybe (|54|) plus 
an integral term with Levy measure as in [5] . There is no general restriction 
but the idea is that a certain generalized symbol function is bounded on 
its domain of definition. More precisely, recall that the ordinary symbol is 
defined by 



a : $7 x il 



C, 



G (X, IU) 



jk 



+ / [exp(iuz) — 1 — iu ■ z1d(z)] v(x, dz) 



for some Q C M. d . Furthermore, the symbol function associated with the 
operator Aq is 



ctq : Q x 



C, 



a (x, iu) = ~ X a °jk( x ) u j u k + X b °j' 



X )IUa 



jk 



+ / [exp(iuz) — 1 — iu ■ z1d(z)] v°(x, dz) 



Now asuume that the Cauchy problem with operator Aq has an explicit 
solution in the form (|55p . First we define the generalized symbol function 
by 



(j g : f2 x il 



a g (x,iu) := a(x,ip (t,u)). 

Here the partial derivative d\ is the partial derivative with respect to the 
first variable, i.e. the time variable. Furthermore and similar as in the case 
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of ordinary symbols we define 

a g i : iR d -> C, 

(T g i(iu) := (Ti(i/j (t,u)), 

where 07 is defined as before. Now there is another interesting effect: using 
the fact that 

d t exp (<j>o(t) + ip (t, iu)x) = A exp (4>o(t) + ip (t, iu)x) (57) 

consider the ansatz p(t,x,u) = exp (4>o(t) + i/jo(t,iu)x) (1 + d k (x, ipo(t, u))t k ) 
and plug it into the Cauchy problem 

d t p = Ap (58) 

Similar computations as before (using the fact that d k has only powers of 
the spatial variable up to order A;) leads to the recursion 

do = 1 

4+i = fpt £i<| e |<fc K 3 ^ ^o) - ff o(aJ, ifo)) d k (59) 
+ kTiJ2o<\e\<k d l (7 ( x ^o)^d e x d k 

Note that if do = a, then d k = for k > 1 as espected. Now the point is 
that we have indeed a bounded symbol, if ip is bounded. Well it makes sense 
to rewrite the latter recursion in the following way 

4 = 1 

4+i = kTi^i<\e\<k Aa ( x ^o)d k + 

(60) 

kh ^o<\e\<k 3f Aa(x, Vo)^<%4+ 

kTl Si<|6|<fc ^M x > ^o)?d x d k , 

where Aa(x,ip) = a(x,ijjo) ~ &o{&) ^o)- Note that in the second sum the 
multiindex starts for multindices e with |e| > 1. For example, we may 
have an unbounded symbol where the unboundedness is caused by the drift 
terms. Then the dmrence symbol is bounded and all derivatives of the 
symbol function are bounded and hence the symbol expansion holds in the 
whole domain of the symbol variable £. Note that this point could not 
be made if we just consider a brute force expansion for generalized symbol 
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functions via the recursion 

4+1 = fc^r ({-dtfoit, u) - xdi^ (t, u))d k 

(61) 

+ T,\e\<k d \ (J ( X ^)^ d l d k), 



where d\ denotes the partial derivative with respect to the time variable 
t. Results analogous to theorem 2.11 and theorem 2.13 (or thorem 2.5. or 
theorem 2.9 in the univariate case) can be written down for either recursion 
(|60|) or (|6ip . As an example one may consider the univariate case where the 
solution of the equation 

- = (a + a 1 x)-^ + (b + b 1 x)- 

with initial data p(0,x, u) = exp(iux) can be easily found via solving the 
corresponding Riccati equations. Not in all cases these solutions correspond 
to an affine process, of course (but we do not need this assumption, if we 
consider the solutions to be just solutions of Cauchy problems which have 
an interpretation to be characteristic functions sometimes). Or one may 
consider stochastic volatility models as in [5] where the explicit solution of 
the corresponding Riccati equations is well known, i.e. Cauchy problems 
determined by the stochastic differential equation 

dX t = (b 10 + b n v)dt + y/vffidW?, 

(62) 

dv = (6 20 - b 2 iv)dt + a^/v{€)dWl 



where W{ and are standard Brownian motions with constant correlation 
p Heston determined the characteristic function to be 

p{t, x, v; u) = exp (<f>(t, u) + ^01 (^ u )v + iux) , (63) 

where 

<f>(t, u) = b 00 iut + ^ { (6 2 i - paiu + d)t-2 In 



(64) 



7T 



1— exp(dt) 
1— gcxp(dt) 



and 



621 — poiu + d 



9=7 : 1, d = sj (paiu - &21) 2 - cr 2 (2b u iu - u 2 ) (65) 

621 — p&iu — d 

In this case the generalized symbol function would be of the form 

a g (x,iu) = a(x,iJjo), where tpo = (ipoi,iu) (66) 
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Note that such an expansion with generalized symbols may be useful also 
from a numerical/computational point of view - sometimes we may have 
faster convergence but this has to be checked from case to case. Furthermore, 
we note that our formulas lead to some results about global solutions of 
Cauchy problems on the flat n-torus for operators with affine coefficients 
which are not covered by the general results of Hormander. Only some of 
that results have an probabilistic interpretation. 

5 Further remarks 

Generalizations of our formulas for the characteristic function to more gen- 
eral Feller processes seems difficult (if not impossible). In can be shown 
that on the flat <i-torus divergence phenomena for expansions considered in 
this paper are quite persistent (cf. [3] for the univariate case). A way out 
may be to consider other holomorhic transforms (therefore a more general 
context is considered in [3]). A refined analysis of the affine counting tri- 
nagles may reveal that for some processe the power series provided in this 
paper converge also on unbounded spatial domains. Implementation and 
numerical analysis of the formulas presented in this paper will be considered 
in a subsequent paper. From the computational point of view this is very 
interesting because we can avoid the difficult numerical solution of general- 
ized Riccati equations. A. Eisenblatter (formerly ZIB, Berlin) pointed out to 
me how the integer recursions considered in this paper may be programmed 
efficiently using rather sophisticated languages like Haskell etc.. Well this 
will be considered in a subsequent paper. 
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